Experimental observables near a nematic quantum critical point 
in the pnictide and cuprate superconductors 
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The newly discovered high temperature superconductor SniFeAs(Oi-xFx) shows a clear nematic 
transition where the square lattice of Fe ions has a rectangular distortion. Similar nematic ordering 
has also been observed in the cuprate superconductors. We provide a detailed theory of experimental 
observables near such a nematic transition: we calculate the scaling of specific heat, local density of 
states (LDOS) and NMR relaxation rate 1/TiT. 



Rapid and important progress has been made in stud- 
ies of Iron-oxypnictides superconductors. Various sam- 
ples with similar FeAs plane and rare earths have been 
synthesized and several different compounds have shown 
superconductivity over 50K [l|, y, y, y, |5|, ll7| , when the 
parent compound is doped by either fluorine or oxygen 
deficiency. Although many experimental facts, including 
the pairing symmetry, are still under debate, all these 
different samples share common features: a tetragonal- 
monoclinic (orthorhombic) lattice distortion and (tt, 0) 
spin density wave (SDW) which commonly exist in the 
undoped samples, and "compete" with superconductiv- 
ity at finite doping. The SDW and lattice distortion 
were first observed in LaFeAs(Oi_xFx) by elastic neu- 
tron scattering and X-ray spectroscopy [6| . Later on this 
phenomenon was confirmed in many other samples with 
La replaced by Ce [7|, Sm 



samples ^ 
{lili^, and 

SrFeaAsa 



17, m and Nd 



also in oxygen free materials BaFe2As2 [S, _j, 
illll], and CaFe2As2 [il,IIl|. In all the samples, at the 
lattice distortion temperature Td, the resistivity shows 
a A shaped anomaly; therefore the A anomaly of resistiv- 
ity can be taken as a measure of the lattice distortion in 
experiments. 

Both lattice distortion and SDW are suppressed under 
doping, but, in general, the lattice distortion occurs at 
higher temperature than the SDW. In SmFeAsOi_xFx, 
the lattice distortion and superconductivity coexist in 
a finite range of doping; the lattice distortion temper- 
ature seems to vanish within the superconducting phase 



171 Il8l | , while the coexistence between SDW and super- 



conductor was never observed. In Ref. [2lJ, |22|, the 
lattice distortion is attributed to anisotropic antiferro- 
magnetic correlation between electrons along x and y di- 
rections, without developing long range SDW. Since this 
order deforms the electron Fermi surface, equivalently, it 
can also be interpreted as electronic nematic order. This 
nematic order has Ising symmetry, therefore the transi- 
tion temperature is controlled by the intralayer spin in- 
teraction, while the long range SDW is controlled by the 
interlayer spin interaction which is much weaker. There- 
fore the nematic transition (lattice distortion) occurs at a 
higher temperature than the SDW in general, and unless 
very close to the critical point, the nematic transition at 



finite temperature should belong to the 2d Ising univer- 
sality class. The distance between the lattice distortion 
temperature and the SDW temperature depends on the 
anisotropy between ab plane and c axis, which can be 
checked by comparing the anisotropy of different sam- 
ples. 

The intimate relation between the structure distortion 
and SDW phase proposed by Ref. 2l|, [221 has gained 
support from recent experiments. It is suggested by de- 
tailed X-ray, neutron and Mossbauer spectroscopy stud- 
ies that both the lattice distortion transition and the 
SDW transition of LaFeAs(Oi_xFx) are second order 
[23|, where the two transitions occur separately. How- 
ever, in AFe2As2 with A = Sr, Eu, Ba, Ca, the struc- 
ture distortion and SDW occur at the same temperature, 
and the structure distortion becomes a first order tran- 
sition 10|, 12, [isl, |lj, |l6| (or a very steep second order 
transition 24!|). These results suggest that the SDW and 
structure distortion are indeed strongly interacting with 
each other, and the structure distortion is probably in- 
duced by magnetism. 

We also note that the importance of nematic ordering 
has also been discussed recently in the context of the 
cuprate superconductors 25|, |20| . Our results below are 
presented in the context of the pnictides, but all of the 
scaling properties of the experimental observables apply 
equally to the cuprates. 

We focus on the zero temperature nematic phase tran- 
sition at finite doping, motivated by the experimental 
suggestion of the existence of structure distortion crit- 
ical point within the superconducting phase of sample 
SmFeAsOi_xFx Hill- By contrast, the SDW phase 



shows no overlap with the superconducting phase in all 
the samples studied so far, therefore we will generally 
ignore it except for noting that the transition from the 
SDW to the nematic order is likely an z = 1 0(3) transi- 
tion, based on the fact that the SDW order wave vector 
is independent of doping [7| , so the low energy particle- 
hole excitations at the SDW wave vector vanishes rapidly 
with small doping and hence make no contribution to the 
damping of the SDW order parameter 21| . The univer- 
sality class of the nematic transition strongly depends on 
the pairing symmetry of the superconducting phase. If it 






FIG. 1: The global phase diagram. The blue, red and green 
curves are phase boundaries of SDW, nematic (structure dis- 
tortion) and superconductivity respectively. 



is an s— wave superconductor without nodes, the transi- 
tion of nematic order will be an ordinary 3D Ising transi- 
tion; while if the superconductor is d— wave, the gapless 
nodal particles may change the universality class of the 
nematic transition. The recent STM [27| and Andreev re- 
flection measurement [28'| suggest that SmFeAsOi-_xFx 
has nodes in the cooper pair, and the spin susceptibil- 
ity measured by Knight shift will tell us whether it is a 
p— wave or d— wave pairing. In our current work we as- 
sume a d— wave pairing. The universal behavior of the 
nematic transition in a d— wave superconductor was first 
studied in Ref. [25]. Recently the same theory was stud- 
ied carefully, and it was shown that in the infrared limit 
there is a special fixed point with logarithmically diverg- 
ing velocity anisotropy of the nodal particles [2f| . In the 
current work we will calculate experimentally relevant 
quantities close to this nematic quantum critical point. 
The global phase diagram is shown in Fig. [TJ 

The low energy Lagrangian describing the nematic or- 
der and nodal particle reads [25j 



FIG. 2: The Feynman diagram used in this work, the dashed 
lines are the propagators of nematic order parameter 0. a, the 
self-energy correction to fermion ^P; b, the self-energy correc- 
tion to 0; c, the vertex correction to fermion bilinear ^tr^vfj; 
d, the vertex correction to fermion bilinear ^^Ta^j- 




FIG. 3: Plot of local density of states with hao/^/o ~ 1/5, 
yAo/vfo ~ 1/10, VAo/vfo = 1/20 from the top to bottom. 
The horizontal axis is in scale of uj/Tc- The curves are y = 
x"-^^ (red), y = x"-^'"'' (blue), y = x°-^'''' (green), y = x 
(black). 



^(p + ^^(p, 



Nf 



-^* = H *la(^r - iVfd^T^ - iWAaj,T^)*la 



a=l 



+ *L(^r - iVfdyT'' - iVAdxT'')^2a, 



L, = Ud.^Y + '-(y^f + U^ + ^^\ 



c 
"2 



4! 



Nf 



i*0 = Ao0^(*Lr"*la + *L^"*2a). 



(1) 



The Nambu fermion ^ is defined in the standard con- 
vention: *ia = (/la, Cafc/gb) and *2a = {ha, <^abflb)- « 

and b are spin indices, /i, /2, fs and f^ are slow fermion 
modes at nodal points {Q,Q), {—Q,Q), {—Q,—Q) and 



{Q, —Q) respectively. If the system develops long range 
order of (/), the four nodal points of the d— wave supercon- 
ductor will be shifted and break the Civ symmetry down 
to C2v due to the couphng Lij-^ [25[. The lagrangian ([1]) 
is not Lorentz invariant because in the real system v^/vf 
is in general not unity. Also, the coupling L^^ breaks the 
Lorentz invariance, since \|/'I't-^\1/ is only one component 
of the space-time current of the Dirac fermion. Therefore 
a realistic scaling procedure is to allow wa and Vf flow 
independently under renormalization group (RG). 

The nematic transition fixed point in Ref. [26| was ob- 
tained by expansion of VA/ivfNj), assuming a small ini- 
tial value VAo/vfo, and in the current situation Nf = 2. 
The RG flow of velocities is obtained by calculating the 
one-loop self-energy in Fig. [^h, with dressed propa- 
gator at order 1/Nf in Fig . ^. The renormalization 



condition is chosen to be keeping the couphng constant 
Aq in L^rf, invariant under RG flow. After the one-loop 
correction, the flow of the self-energy and velocities reads 



din A 

dvf 
dlnA 

duA 
dlnA 

d{vA/vf) 



= Ci{-iu) + C2VfkxT'' +C3VAkyT'', 

= iC,-C2)vf, 

= {Ci-C3)VA, 

= {C2 - C3){VA/Vf). 



(2) 



dlnA 

A is the momentum cut-off. Ci , C2 and C'3 arc functions 
of v^/vf and Nf, their detailed forms are given in the 
appendix. Using these RG equations, we are ready to 
calculate the scaling of the local density of states (LDOS) 
accessible by scanning tunneling microscope (STM): 



p{u) 



QjK'r GiKi 



^Tr{ImGi,„t(u/fc:r, WA^y, w) 



(27r)2 

linG2,ret{vfky,VAkx-,t^)] 



VAVf 

dkLdk' 



——^lmGi^ret[kx,ky,Uj) 



{2nf 



'-luiG2.xet{k'y,k'^,Uj)}. 



(3) 



Gi and G2 are retarded single particle propagator for 
\l/i and ^2- The RG equations in Eq. [2] are calculated 
by rescaling momentum cut-off. Since va/vj- is much 
smaller than 1 and flows to zero under RG, for frequency 
w, the corresponding momentum scale is p = ^/vf. 
Therefore dhiLo/dhip = 1 + dhivf /dhip. Now the scal- 
ing equation for p{uj) reads: 



dlnp 
dhibj 



din p 



dlnp 



dlnpx ^f^ 

^ amp 



'^inAxa-^^; 



a-ci)- 



din 



v"/"A. 



din A 



1- 



d In 1 



din A 



1 — G3 — G2 -l- Ci 

1 + C2-CI ■ 



(4) 



The ultraviolet cut-off of the theory is taken to be the 
transition temperature Tc at the critical doping of ne- 
matic transition. Although va will be renormalized to 
be zero in the infrared limit, the expansion of G2 and G3 
with small va/vj given in Ref. [261] shows that va ap- 
proaches zero slowly with energy scale, therefore for the 
experimentally relevant energy scale, one cannot naively 
take the flxed point value of va- Instead, we have to 
integrate Eq. |4] numerically from the ultraviolet cut-off, 
and the result at certain frequency oj depends on the ini- 
tial value of vao/vfo- The results of p{uj) with initial 




FIG. 4: Plot of nodal fermion contribution to specific heat 
with VAo/vfo = 1/5, VAo/vfo = 1/10, VAo/vfo = 1/20 from 
the top to bottom. The horizontal axis is in scale of T/T^. 
The curves are y = x^-^^ (red), y = 2;^''^ (blue), y = a:^'^"^ 
(green), y = x^ (black). 



value vao/vfo — 1/5, 1/10, 1/20 are plotted in Fig. [31 
for frequency between e~^ < oj/Tc < e^^ . unlike ordi- 
nary d— wave superconductor with LDOS p{i^) ~ w, in 
all three plots the LDOS scales with frequency as 



p{uj) 



a < 1. 



(5) 



This simple power law relation fits well for the frequency 
range e^^ < uo/Tc < e~^, and it can be checked by STM 
technique on samples in the quantum critical regime. 

The fluctuation of nematic order parameter certainly 
affects the thermal dynamical quantities. The free energy 
J^ — TlnZ/V has scaling dimension d + z, the singular 
part of the free energy can be written as 



^F^iUxCyY 



(6) 



Now temperature is taken to be the infrared cut-off: i^^- ^ 
1/T, and the spatial correlation length can be estimated 
to be ^x '^ VxS,Ti S.y ~ VyS^r- The anisotropic velocity of 
^i leads to the following contribution to the free energy: 



J'-i! 



dlnG* 
dlnT 



1 



VfVA 



rr>6 



din I 



"/"A. 



dlnT 



1 + G2 - Gi 



- 2Gi — G2 — G3 
1 + G2 - Gi 



(7) 



The (f) contribution to the free energy and speciflc heat 
can be estimated in the same manner, although there is 
no velocity anisotropy for the (j) fleld. The velocity of (fi 
in the large Nf limit can be evaluated by calculating the 
one loop correction to the self-energy of (j). In the case 
of small VA/vf, the velocity of can be taken to be Vf 



isotropically: 



d In C^ 
d\nT 






2 + 2Ci - 2C2 
I + C2-C1 ■ 



(8) 



The solutions of Eq. [7]with different v^q/vfq are plotted 
in Fig. [4] The equations are solved for the experimentally 
relevant temperature range e"** < T/Tc < e~^. In gen- 
eral, the nodal ferniions contribute more to the specific 
heat compared with field, because va scales stronger 
with temperature compared with Vf. For e~^ < T/T^ < 
e~^ , the scaling of specific heat is 



C ^T^, /3 < 2, 



(9) 



which distinguishes the current situation from the ordi- 
nary d— wave superconductor with nodes. 

Another way to probe the density of states is the NMR 
relaxation rate l/Ti, which is related to the following 
Green function: 



F{oj) 



dq^dqy-x"{qx,qy,uj), (10) 



in the limit of w — > 0. The scaling of l/(TiT) with tem- 
perature is the same as the scaling of F{lo) with w, as 
T and uj can both serve as infrared cut-off of the theory. 
The momentum integrated susceptibility should involve 
spin density at various "slow" momenta. At the low en- 
ergy theory of the nodal particles, following fermion bi- 
linears are low energy spin density modes that have uni- 
versal scalings: 

q = (0,0), *lcr'^*i-|-*t^'^^2; 

qi2A = (2g,0), ^t^'^*!; 

-qi2A = (-2Q,0), ^la-^^a; 

qi2B = (0,2Q), **tW^*2; 

-qi2B = (0,-2g), *^r%V^**, 

qn = (2g,2Q), ^^rV^-^^i; 

-gn = (-2Q,-2Q), vptrW**, 

922 = (-2Q,2g), **2tW*2; 

-922 = (2g,-2g), ^Irya-ay^; 



(11) 



ct" are three spin Pauli matrices. To evaluate F{lu) we 
need to calculate the correlation of all the fermion bilin- 
ears above. The susceptibility x gains fermion self-energy 
correction as in Fig. [2ji as well as vertex corrections Fig 
H: and Fig. ^ for vertices *|Ta*j and **rA*j respec- 
tively. For a general fermion bilinear with flavor matrix 




FIG. 5: Plot of l/{TiT) with T from contribution of spin 
density at q = (0,0), with VAo/vfo = 1/5, VAo/vfo = 1/10, 
VAo/vfo = 1/20 from the top to bottom. The horizontal axis 
is in scale of T/Tc- The curves are y — x^'"^"^ (red), y — a;^'*'^ 
(blue), y = x^'^'^ (green), y — x^ (black). 



Ta, the vertex correction RG equation is conventionally 
written as 



dTA 
dlnA 



CaT 



A^A- 



fl2) 



Ca is a function of VA/vf. The details of calculations are 
given in the appendix, the results are 



Co 



Cr 



"Ol, Ljj-j: 



Ci 



^Cs 



C3 — Ci — C2, C'r^ 



-C2 



Ci2A = Ci2B = -0.3486 

Cii = C22 = ^C'rv- 



Nf 



(13) 



Co is the vertex correction for fermion bilinear '^j'^i 
and ^Jcr^^i, because spin is a good quantum number 
spin Pauli matrices do not change the vertex correc- 
tion. Ci2yi and Ci2B are vertex corrections to \I'J^ct°^2 
and \l/^T^cr^cr"\I'2, Cn and C22 are vertex corrections to 
^'iT^o-^ct"*! and *^r''cr''a°*2 respectively. 

After taking into account of both self-energy and ver- 
tex corrections, the scaling equation of the relaxation rate 
reads: 



d\nF{uj 
dlnuj 



2 - 2Ci - 2Ca + 2- 



dln 



^"/"A/ 



d In c 



I + C2-C1 



2 - 2Ca - 2C2 - 2C3 + 2Ci 
I + C2-C1 ■ 



(14) 



Equation (fT4| leads to the following scaling equation for 
l/{TiT) with temperature: 



din l/{TiT) 2-2Ca- 2C2 - 2C3 + 2Ci 



dlnT 



1 + C2 - Ci 



(15) 




FIG. 6; Plot of l/{TiT) with T from contribution of spin 
density at q — (2(5,0) and q— (0, 2Q), with uao/w/o = 1/5, 
VAo/vfo = 1/10, VAo/vfo = 1/20 from the top to bottom. 
The horizontal axis is in scale of T/Tc- The curves are y — 
x^-^^ (red), y = x^'^ (blue), y = x^'^'^ (green), y = x'^ (black). 




FIG. 7; Plot of l/{TiT) with T from contribution of spin 
density at q= {2Q,2Q) and q= {2Q, — 2Q), with VAo/vfo = 
1/5, VAo/vfo — 1/10, VAo/vfo ~ 1/20 from the bottom to the 
top. The horizontal axis is in scale of T/T^. The curves are 
y = a;2"» (red), y = x'^-°^ (blue), y = x^"^^ (green), y ^ x'^ 
(black). 



The contribution from different fermion bilinear compo- 
nents listed in Eq. [Tl] should be solved individually, and 
the solutions of each component are plotted in Fig. [51 
[S] and [T] One can see that at small frequency the spin 
density at (0, 0) makes the most substantial contribution, 
which makes the scaling of l/(TiT) differ from the ordi- 
nary d— wave superconductor: 



TiT ' ' 



(16) 



The superconductivity can be fully suppressed by 
strong enough magnetic field. Since the nematic quan- 
tum critical point is most likely in the underdoped regime 
with low Tc, it is possible to apply strong enough inplane 
magnetic field H > Hc2 in experiments. With supercon- 



ductivity fully suppressed by inplane magnetic field, the 
universality class of the nematic transition becomes the 
z — 3 theory described in Ref. 2]J. The z — 3 quan- 



tum critical point leads to a large number of low energy 
excitations, which will contribute to thermal dynamics 
and transport. The standard mean field theory leads to 
following results at low temperature 2l| : 



C-T 



2/3 



y4/3^ 



These results will finally crossover to d 
enough to the critical point: 

C - T, p - T^/^. 



(17) 
: 3 scalings close 



(18) 



Similar z — 3 nematic quantum critical point was dis- 
cussed in Ref. [29] . 

In a summary, in this work we computed the scaling 
of physical quantities close to a nematic quantum critical 
point in a d— wave superconducting phase of the newly 
discovered material SmFeAsOi_xFx, motivated by recent 
experiments on the polycrystal sample. For the exper- 
imentally relevant energy range, the scaling of LDOS, 
specific heat, NMR relaxation rate all deviates from an 
ordinary d— wave superconductor. These results also ap- 
ply, essentially unchanged, to the cuprate superconduc- 
tors |25l . |26]. This research is supported by the NSF 
under grant DMR-0537077. 



APPENDIX 

In this section we shall calculate the vertex corrections 
to the fermion bilinears listed in Eq. [TT] For a gen- 
eral vertex "^JTa'^'j, the vertex correction is calculated 
according to Feynman diagram in Fig. [2t: 

^ ^ 1 /" (P'V dflr „„ ,„ 

T^ = Ta + — I 77TTT2 — [^ G,(f},p) 



Nf J (27r)2 2tt 
xTAG,{n,p)T' 



1 



T2in,p) 



(19) 



For a general vertex ^*Ta$j, the vertex correction is 
calculated according to Feynman diagram in Fig. ^d- 



rA-TA + ^ f ^^i-l)[r^Gdn,p) 



Nf J (27r)2 27r 

xTAG,{n,p)T-] 



1 



r2(r!,p)- 



(20) 



r2(r2,p) is the self-energy of (j) from integrating out 
fermions: 



T2iuj,p) = Il2{kx,ky,uj)+Il2{ky,kx,uj) 

1 ^2 _|_ „2^2 



U2{kj;,ky,Uj) 



IGvfVA (w^ + vjkl + wi/c^)!/ 



T7J-(21) 



112 is reminiscent of the vacuum polarization of the 2+ Id 
QED, with a gauge invariant form Il^i,{p) '^ p{S^i, — 
PfiPv/p^), and Ha = Tl^xip)- 

The RG equation is defined as the change of parame- 
ters after rescahng the cut-off. The cut-off is introduced 
in a smooth function /C(p^/A^), with A^(0) — 1 and falls 
rapidly when y — > 1. Now following the notation in Ref. 
[26|, we change the momentum space integral to cylin- 
drical coordinates p^ = yK{v fX,cos9 ,s,m.9), and the RG 
equations in Eq. [2] and Eq. [12] are obtained when the 
cut-off is reduced from A to A — dA, with 



Ci — —Co — 



2(fAM) 



dx I de 

oo Jo 



(cc^ — cos^ 9 - {vA/vfY si 



sm 



C, 



(x^ + cos^^ -h (i'a/^'/)^ sin i 

-c.. ^'M^ r dx r d9 



;Q{x,9)- 



(— a;^ + cos^ 9 - (va/vj) 



^sin^ 



(x^ -I- cos^ 9 + [vA/vfY sin 



2fl^2 



■Q{x,9)- 



-2(-Mrd.rd. 

/ J -oo Jo 

,2 I „„„2 n (.. I.. \2 „:„2 



O3 0^-3= 

(a;^ -I- cos^ 9 — [vaIvjY sin 
(x2 + cos2 9 + (wa/w/)^ sin^ 6*) 2 



■a(x,0); 



a 



-Cii — —Ci 



2{vA/vf) 



V/ 



dx I d9 
Jo 



{x''+cosH+ivA/vf?smH) 
{x^+cos'^9+{vA/vfysm^9y^ ' ' 

= C3 — Ci — C2; 



Ci2A — Ci2B — 

i-x') X 



2{vA/vf) '■+°° '■2" 



1 



dx / d6i 





g-'{x, 



x^ + cos2 9 + [vA/vfY sin^ 
x'^ + {vA/vfY cos'^9 + siT? 9 



x^ + sin^ 6* 

J x"^ + {vA/vfY cos^ + sin^ 
x^ -f' cos^ 9 



va/vjY sit? 9 A- cos^ 9 



(22) 



It was noted in Ref. [26| that function C3 has a rather 
singular form in the small va/v/ limit: 



C3 = ^[ln(./M)- 0.9601]^ -I- 



Plugging this function into Eq. [2l one can see that 
va/vj approaches zero a little faster than the ordi- 
nary marginally irrelevant operators, but for the exper- 
imentally relevant energy scale, va/v/ still flows slowly. 
Therefore all the plots in our paper, though integrated 
from a complicated equation, can be fit with a simple 
power law. 



[1] 
[2] 

[s; 

[4; 
[5; 

[6; 

[7 
[8 

[lo; 

[11 

[12 
[13 

[14; 

[15 
[16 

[17 



(23) 



Linjun Li, Yuke Li, Zhi Ren, Xiao Lin, Yongkang Luo, 

Zengwei Zhu, Mi He, Xiangfan Xu, Guanghan Cao, 

Zhu'an Xu, arXiv:0806.1675 (2008). 

Zhi-An Ren, Jie Yang, Wei Lu, Wei Yi, Guang-Can 

Che, Xiao-Li Dong, Li-Ling Sun, Zhong-Xian Zhao, 

arXiv:0803.4283 (2008). 

Zhi-An Ren, Guang-Can Che, Xiao-Li Dong, Jie Yang, 

Wei Lu, Wei Yi, Xiao-Li Shen, Zheng-Cai Li, Li-Ling 

Sun, Fang Zhou, Zhong-Xian Zhao, EuroPhys. Lett. 83, 

17002 (2008). 

X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, D. F. 

Fang, Nature 453, 761 (2008). 

Linjun Li, Yuke Li, Zhi Ren, Xiao Lin, Yongkang Luo, 

Zengwei Zhu, Mi He, Xiangfan Xu, Guanghan Cao, 

Zhu'an Xu, larXiv:0806.1675l f2008'). 

Clarina de la Cruz, Q. Huang, J. W. Lynn, Jiying Li, W. 

Ratcliff H, J. L. Zarestky, H. A. Mook, G. F. Chen, J. 

L. Luo, N. L. Wang, Pengcheng Dai, Nature 453, 899 

(2008). 

Jun Zhao, Q. Huang, Clarina de la Cruz, Shiliang Li, J. 

W. Lynn, Y. Chen, M. A. Green, G. F. C hen, G . Li, Z. Li, 

J. L. Luo, N. L. Wang, Pengcheng Dai, larXrv :0806.258g1 

(2008). 

Marianne Rotter, Marcus Tegel, Dirk Johrendt, 

arXiv:0805.4630 (2008). 

N. Ni, S. L. Bud'ko, A. Kreyssig, S. Nandi, G. E. Rustan, 

A. L Goldman, S. Gupta, J. D. Corbett, A. Kracher, P. 

C. Canfield, arXiv:0806.1874 (2008). 

Q. Huang, Y. Qiu, Wei Bao, J.W. Lynn, M.A. Green, 

Y. Chen, T. Wu, G. Wu, X.H. Chen, arXiv:0806.2776l 

(2008). 

Marianne Rotter, Marcus Tegel, Inga Schellenberg, 

Wilfried Hermes, Rainer P02ttgen, Dirk Johrendt, 

larXiv:0805.4 021 (2008). 

C. Krellner, N. Caroca-Canales, A. Jesche, H. Rosner, A. 

Ormeci, C. Geibel, arXiv:0806.1043 (2008). 

J.-Q. Yan, A. Kreyssig, S. Nandi, N. Ni, S. L. Bud'ko, 

A. Kracher, R. J. McQueeney, R. W. McCallum, T. A. 

Lograsso, A. L Goldman, P. C. Canfield, larXiv:0806.27lTI 

(2008). 

Jun Zhao, W. Ratcliff H, J. W. Lynn, G. F. Chen, 

J. L. Luo, N. L. Wang, Jiangping Hu, Pengcheng Dai, 

irXiv:0807.1077 (2008). 

Milton S. Torikachvili, Sergey L. Bud'ko, Ni Ni, and Paul 

C. Canfield, arXiv:0807.0616 (2008). 

A.L Goldman, D.N. Argyriou, B. Ouladdiaf, T. Chat- 

terji, A. Kreyssig, S. Nandi, N. Ni, S. L. Bud'ko, P.C. 

Canfield, R. J. McQueeney arX iy:0807, l525 (2008). 

R. H. Liu, G. Wu, T. Wu, D. F. Fang, H. Chen, S. Y. Li, 

K. Liu, Y. L. Xie, X. F. Wa ng, R. L. Yang, L. Ding, C. 

He, D. L. Feng, X. H. Chen, larXi^:0804.2105 (2008). 



[18] Serena Margadonna, Yasuhiro Takabayashi, Martin T. 

McDonald, Michela Brunelli, G. Wu, R. H. Liu, X. H. [24] 
Chen, Kosmas Prassides, arXiv:0806.3962 (2008). 

[19] Y. Qiu, Wei Bao, Q. Huang, J.W. Lynn, T. Yildirim, J. 

Simmons, Y.C. Gasparovic, J. Li, M. Green, T. Wu, G. [25] 

Wu, X.H. Chen, arXiv:0806.2195 (2008). 

[20] Ying Chen, J. W. Lynn, J. Li, G. Li, G. F. Chen, J. L. 

Luo, N. L. Wang, Pengcheng Dai, C. dela Cruz, H. A. [26] 
Mook, arXiv:0807.0662 (2008). [27] 

[21] Cenke Xu, Markus MueUer, Subir Sachdev, Phys. Rev. 
B 78, 020501 (R) (2008). 

[22] Chen Fang, Hong Yao, Wei-Feng Tsai, JiangPing Hu, [28] 
Steven A. Kivelson, Phys. Rev. B 77 224509 (2008). 

[23] M. A. McGuire, A. D. Christianson, A. S. Sefat, B. C. [29] 
Sales, M. D. Lumsden, R. Jin, E. A. Payzant, D. Man- 
drus, Y. Luan, V. Keppens, V. Varadarajan, J. W. Brill, 
R. P. Hermann, M. T. Sougrati, F. Grandjean, G. J. 



Long, arXiv:0806.3878 (2008). 

Marcus Tegel, Marianne Rotter, Veronika Weiss, Falko 

M. Schappacher, Rainer Poettgen, Dirk Johrendt, 

arXiv:0806.4782 (2008). 

Eun-Ah Kim, Michael J. Lawler, Paul Oreto, Subir 

Sachdev, Eduardo Fradkin, Steven A. Kivelson, Phys. 

Rev. B 77, 184514 (2008). 

Yejin Huh, Subir Sachdev, arXiv:0806.0002 (2008). 

Oded Millo, Itay Asuhn, Ofer Yuli, Israel Felner, Zhi-An 

Ren, Xiao-Li Shen, Guang-Can Che, Zhong-Xian Zhao, 

arXiv:0807.0359 (2008). 

Yonglei Wang, Lei Shan, Lei Fang, Peng Cheng, Cong 

Ren, Hai-Hu Wen, arXiv:0806.1986 (2008). 

Vadim Oganesyan, Steven Kivelson, Eduardo Fradkin, 

Phys. Rev. B 64, 195109 (2001). 



